Theory of Single Photon Control from a Two Level System Source 
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Generation of a single photon or a pair of photons from a single emitter is important for quantum 
information applications. Using the generating function formalism we investigate the theory of a 
few photons on demand for the square laser pulse and the rapid adiabatic following method. Exact 
theory and numerical solutions are used to design control fields for a two level emitter, which yield 
an optimal single or two photon source, under the constrains of finite laser field strength and finite 
interaction time. Comparison to experiments of Brunei et al, shows that the experiments were made 
close to optimal conditions. 
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PACS numbers: 42.50.-p, 33.80.-b,32.50.+d 



The generation of a single photon on demand, is mo- 
tivated by fundamnctal interest in quantum properies of 
light as well as possible applications like quantum 
computation, cryptology, and communication @>E| ■ Sin- 
le photon experiments using two level rubidium atoms 
, single organic molecules H, IE 0j an d single quan- 
tum dots §], are nowadays considered in detail due to 
these applications. Besides single photon sources, sources 
of pairs of photons are also important, due to useful 
properties of indistiguishable photons and entanglement 
P3- HE 113 • A generic problem in single photon sources, 
are fluctuations in the number of photons emitted, and 
in the time of emission, due to the quantum nature of 
the emission process. 

Let Pisf(t) be the probability of N emission events 
in the time interval (0,i). The generating function 
2y(s,t) = J2n=o sN Pn {t) contains information on pho- 
ton statistics necessary for the determination of Pjv(i) 
fl2l| . For a two level atom or molecule, within the rotat- 
ing wave approaximation, the generalized optical Bloch 
equations which yield the generating function y(s, t) are 

U(s,t) = -\U (s,t) + S(t)V (s,t) 

V (s, t) = -5{t)U (s, t) - § V (s, t) - n(t)w (s, t) 

W (s, t) = Q(t)V {s, t) - 1 (1 + a) W (s, t) - E (1 + s) y (s, 



experiments. The first is the rapid adiabatic following 
method 6], in this method the single molecule is inter- 
acting with a continuous wave (cw) laser field, and at 
the same time the absorption frequency of the molecule 
in modulated periodically using a slow radio frequency 
Stark field. In this case the Rabi frequency fl(t) = O 
is time independent, and S(t) — -^p- cos (fRFt) @. The 
time dependent detuning 5(t) is designed to bring the 
molecule in and out of resonance with the cw laser in 
such a way that a few photons are emitted per crossing 
(see details below). In a second very common approach 
the emitter interacts with a sequence of laser pulses. We 
model such behavior by the text book example of a square 
pulse where the Rabi frequency is 



0, < t < T 
T < t, 



(2) 



y(s,t) 



(l-s)W(s,t)-^(l-s)y(s,t), 



(1) 

where T is the spontaneous emission rate, and £l(t),8(t) 
are time dependent Rabi frequency and detuning. One 
goal of this paper is to use Eq. Q and find optimal con- 
trol fields d(t) and fl(t) which under physical constrains 
e.g. not too strong laser fields, yield a few photons. Inter- 
action of the molecule with its environment, for example 
effects of spectral diffusion 0, 0, and other sources 
of stochasticity 18] like triplet blinking, are not consid- 
ered here. Still as we show later the generating function 
formalism is in excellent agreement with cryogenic tem- 
perature single molecule experiments. 

Two types of control approaches are used so far in 



and the detuning 5(t) = 0. The specific goal in appli- 
cations varies, some experiments are interested in single 
photon sources, while others consider two photon sources. 

The goal of theory of single photon control is to find 
the probability of emission of N = 0, 1, 2 photons. And 
then to design the optimal external fields 8(t),Q(t) to 
maximize the probability of single photon emission or 
the probability of the emission of a pair photons. Since 
very little is known on this important control problem, 
i)we consider here simple cases. In particular, the opti- 
mization of Pi and P2 under the constrain of finite time 
and strength of the excitation laser field is investigated. 

We first investigate the square pulse and later consider 
the rapid adiabatic following method. Only the limit of 
long measurement time t > > T is considered, for the sake 
of space. The details on the calculation of Pat will be pub- 
lished elsewhere, this calculation becomes cumbersome 
already for small N. We use dimensionless units with 
r = 1, the notation y = vT — 4£1 2 and x = \/\ — 16£! 2 , 
and find 



Pn 



2y 2 



J 2 ) cosh(^) 



y-l + 2ysinh(^) 



(3) 



2 



P 



[i-v 2 )e 

8y 5 



y (-8 + T + 4y 2 + Ty 2 ) cosh(^) + 2y (4 + T - 2y 2 - T y 2 ) + 2 (-3 + (1 + T) y 2 ) sinh(^) 

(4) 



Pi 



^ J]/ | [T(T + 8)y 4 + (T 2 ~ 28T - 32) y 2 + 96] cosh 



4 [T (T + 4) y 4 - (T 2 + 8P - 8) y 2 - 24] + 2y [(T 2 - 24) y 2 — T (y 2 + 9) + 60] sinh ( ^ ) [> . ( 5 ) 



(N) 



a 2 



(i + 2fi 2 y 



T (l + 2ft 2 



2ft 2 



-3T 

e 4 



T 



2 - 2ft 2 ) cosh a/I - 16ft 2 - + (2 - 140 



sinh (VI - 16ft 2 f ) 



VI - 16ft 2 



(0) 



(N(N-1)) = 



(l-x 2 ) 



2\2 



_*t 12 - (16 + 3T) x - 2 (6 + T) x 2 + Tx 3 
4 {-3 + xfx 3 



-12 - (16 + 3T)x + 2(6 + T)x 2 + Tx 3 [T 2 (-9 - 



(3 + a;) 4 x 3 



r 



„2\2 



32(63 + 5a; 2 ) + 2T(-351 + 30x 2 



(-9 + a; 2 ) 4 



ill 



Using the moments Eqs. 0) we get an exact expres- 
sion for Mandel's parameter Q = ((N 2 ) - (N) 2 )/(N) - 1, 
which classifies deviations from Poissonian behavior. We 
now investigate limiting cases, which explain the rich 
Physical behaviors of Eqs. !0-[Z|l- 
In the strong field limit ft >> 1 



Po 



e 2 cos 2 



ftT\ 



64 



[4 + 2T - (4 + T) cos (SIT)] 



■T [(8 + T) cos(ftT) + AT + 16] 



(8) 



Hence for fixed T, Pjv shows a bounded oscillatory behav- 
ior, related to the Rabi oscillations of the excited state 
population. Note that Pi and P2 are bounded from above 
according to 



Pi < 



P2 < 



3T)e~ 



TT 



(>4 



(5T + 24) . 



(9) 



While for Pi we may set T — * and get Pi = 1 i.e. a ir 
pulse soon to be discussed, the behavior of P2 is more in- 
teresting and we find that we cannot reach the ideal limit 
of emitting two photons with probability one. Instead we 
find that for a 2ir pulse SIT = 2Tm, with n a positive in- 
teger, the maximum of P2 is p^ ax ~ 0.41 which is found 
when T = 2(V61 — l)/5. Thus unlike the maximum of 
Pi which is found for a very short pulse, the maximum of 
P2 is found for a particular finite interaction time, which 
is needed for the production of a pair of photons. This 
maximum is for a fixed large Rabi frequency. The global 
maximum of P2 is found at ft ~ 1.25, T ~ 4.86 and then 
pMax _ namely the global maximum of P2 is found 
for intermediate Rabi frequencies. 

In the mathematical limit of strong and short pulses, 
the product SIT remaining fixed, we have 



lim Pi 

n^oo,T->0 



ftT 



(10) 



Po = 1 — Pi, namely in this limit either one or zero 
photons are emitted per pulse. The well known 7r pulse 
SIT = nir yields the global maximum of Pi which is unity. 

In experiments fields are never of infinite strength and 
interaction time is not zero and then exact results are use- 
ful. A transcendental equation for the value of T which 
yields the extremum of Pi is found using Eq. JHJ 

- 2P + (2 + T) cos(ftT) + 2(4 + T)ft sin(ftT) = 0. (11) 



3 



Since Pi is oscillating in this limit of strong fields, there 
exists an infinite number of solutions of this equation. 
For large fl the maximum of Pi , p^ax is the first root of 
this equation. For — > oo, Eq. gives fism(fiT) = 
hence the value T ~ mr/Q yields the maximum of Pi and 



pMax ^ e -7r/(2n) | ! + _ 



3tt 
8fi 



(12) 



This equation gives the corrections to the ideal 7r pulse 
limit of fi -> oo and T -> which yields P^ ax = 1. 

The intermediate Rabi frequency marks the transition 
between over damped to under-damped behaviors. For 
fi = 1 /2 we have 



Pi = 



T 2 (480 + 160T + 20T 2 + T 3 ) _ z 



7680 



(13) 



We find Pf ax ~ 0.56 when T ~ 6.75. 

The weak field limit where — > and T — > oo in such 
a way that Sl 2 T is called the semi-classical since 



lim Pn 

n^0,T-*oo 



(n 2 T) 



N 



(14) 



and photon statistics is Poissonian. The semi clas- 
sical maximum of Pi is found when fi 2 T = 1 and 
limn-»o,T-»oo Pj Max = e _1 . While for Pi the maximum 
is found when Q 2 T = 2 and then P 2 Max = 2e~ 2 . 

In Fig. ^we show the maximum of Pi obtained using 
Eq. J2J. The Fig. shows for each Rabi frequency the 
value of the maximum of Pi and the duration of pulse 
interval on which this maximum is found. The general 
trend is that for strong excitation fields p^ax _ ^ f or 
a 7r pulse, and then pMax j g decreasing as we decrease 
17, and at the same time the value of T on which the 
maximum is found is increasing. A transition between 
a strong field and semi-classical limit is easy to observe. 
Fig. [3 shows the maximum of P2 . The three regimes 
of strong fields, intermediate, and semi-classical limit are 
clearly observed. The observation that a peak in P2 is 
found for intermediate Rabi frequency, makes this case 
more interesting, and similar behaviors are expected also 
for other types of pulse shapes. 

For two reasons it is interesting to investigate the prob- 
ability of emission of two photons using a short pulse. 
First since ideal single photon emission demands P2 = 
and secondly since inter-arrival times of a pair of photons 
generated from a short pulse are expected to be shorter 
if compared with inter-arrival times of a pair of photons 
generated by a longer pulse, and hence potentially useful 
for sources of indistinguishable photon pairs. In the limit 
T^Owc find 



480 



(15) 



The T 5 behavior in Eq. (|15|l can be explained, by not- 
ing that the field dependence of P2 is expected, namely 
Pq oc f2°, Pi oc il 2 and P 2 oc SI 4 . Physically we obviously 




100 



Q/r 



FIG. 1: The maximum of the probability of emission of a 
single photon in a square pulse obtained from Eq. |gj with 
r = 20MHz. The circles show the strong field approxima- 
tion Eq. (II 12 P - The star gives the asymptotic semi-classical 



weak Rabi frequency behavior limn^o.T^oo Pi ax = 
dark dot is for the intermediate Rabi frequency Q 
predicted based on Eq. (I13H . 
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FIG. 2: The maximum of P2. For weak fields the semi- 



classical maximum of Pc 



'2 is 



2e 



-2 

>Max 



0.27 while in 

the strong field limit the maximum is P2 Vlax — 0.41 when 
TT = 2(v / 61 - l)/5, as predicted in the text. The global 
maximum (i.e. the dark dot) is found for an intermediate 
Rabi frequency. 



expect that P2 oc Vt i T m and m is a positive integer. As 
mentioned in Eq. 110H for a short and strong pulse we 
have P2 = and hence me must have m > 4, otherwise 
the product f2 4 T m does not approach zero in the limit 
ft — > 00 and T — * with flT remaining fixed. And indeed 
we find m = 5 in Eq. (|15f) . 

Mandel's Q parameter characterizing deviation from 
Poissonian photon statistics, is analyzed using Eqs. 10 
in the limits considered so far for the probabilities P/v ■ 



4 



In the limit of long interaction times 



lim Q = - 

T^oo (1 



6ft 2 



2ft 2 ) 2 ' 



(16) 



This result was obtained by Mandel 0] for an atomic 
transition interacting with a cw laser. As expected if 
the interaction time T is long, namely many photons 
are emitted, the familiar resonance fluorescence sub- 
Poissonian behavior is found. 

New behavior is found for finite interaction time T. 
For strong laser fields, ft — > oo we find using Eqs. (U30) 



(17) 



When the square pulse is short and strong, however ftT 
remains finite, we find an interesting discontinuous be- 
havior 





1 - cos (fiT)e-^ 


9 

ST 

+ 3Tcos(ftT)e- — 


2 


1 + T-cos (OT)e-f 





lim Q 

n->oo,T^0 



-sin 2 (^) ftT ^ 27m 
6/7 ftT = 27m, 



(18) 



where n is a positive integer. Of-course for finite though 
large (small) values of ft (T), Q is a continuous func- 
tion. Though for large ft, Q exhibits: (i) very sharp, 
cusp resonances on ftT = 2im which become narrower 
as the mathematical limit of ft — > oo and T — ► is ap- 
proached and (m) the value of Q on the peaks is roughly 
6/7. To understand the origin of the behavior predicted 
in Eq. i fTSjl note that according to Eq. (fTUf we either 
obtain zero photons or one photon in this limit. Hence 
we expect naively Q = —P± = — sin 2 (ftT/2). However 
this approach cannot be used when Pi ~ 0, namely when 
ftT ~ 27m, when (N) ~ and the fluctuations become 
non-trivial. Instead one must consider the limit with care 
and include the 1/ft corrections, which yield the Q ~ 6/7 
law in Eq. JTHJ). 

We now consider the rapid adiabatic following method. 
In Table 1 we show photon statistics obtained using the 
generating function formalism Eq. Q and experiments 
of the group of Orrit (numbers in the brackets in Table 
1). A good agreement is found, indicating that the simple 
two level system approximation of the molecule, is an 
excellent approximation, a conclusion reached previously 
in The question remains what is the ideal choice of 
control parameters for the generation of a single photon? 
And how does this method compare with the square pulse 
approach? 

From the experimental data in Table 1 we see that 
maximum of Pi for ft = 3.2T is Pi = 0.68 found for 



Arf = 88r. We found numerically (data not shown) 
that for ft = 3.2r, i/rf = 3MHz the maximum of Pi is 
Pi ~ 0.69 for Arf — 100r. We see that the experiments 
of the group of Orrit are conducted very close to optimal 
condition for a particular Rabi frequency, in the sense 
that Pi is close to its maximum. Note that in principle 
for very strong Rabi frequency and very large Arf we 
may obtain effectively a 7r pulse, i.e. a single photon 
with probability one. 



^RF Vrf 



ft 



Pa 



Pi 



Po 



50r 3 MHz 3.2T 0.01(0.02) 0.53(0.56) 0.32(0.31) 
88r 3 MHz 3.2T 0.09(0.11) 0.68(0.68) 0.20(0.18) 
130r 3 MHz 3.2T 0.24(0.22) 0.66(0.66) 0.10(0.10) 
160r 3 MHz 5. Or 0.04 0.72(0.74) 0.21 
Table 1: Rapid adiabatic following method, theory ver- 
sus experiment, T = 20MHz. 

A more detailed investigation of the rapid adiabatic 
method will be published elsewhere. We mention that 
the maximum of Pi for the square pulse is slightly larger 
compared with the rapid adiabatic method. More specif- 
ically we searched for the maximum of Pi in the range 
< ft < 20r, < A RF < 600r and u KF = 3MHz to be 
reasonably close to experimental situations (see Table 1), 
and found p^ Lax - ~ 0.82. This maximum should be com- 
pared with with pMa* ~ 1 in Fig. Q]for the square pulse 
and for a similar range of Rabi frequency. We found that 
the rapid adiabatic following method has one possible ad- 
vantage over the square pulse, and that is the absence of 
fast Rabi oscillations in Pi, in the limit of strong fields. If 
Pi is oscillating the generation of a single photon is less 
stable, and in this sense the rapid adiabatic following 
method is better. The absence of Rabi oscillations is due 
to the smooth time dependent detuning which averages 
out the Rabi oscillations. 

Finally, while our two level system approximation 
works very well for certain low temperature single 
molecules |6j and is expected to work well for simple 
atoms Q, extensions of this work for multi level emitters, 
where even the 7r pulse condition for single photon emis- 
sion is not valid Q should be investigated theoretically, 
though then analytical theory is cumbersome. While our 
work provides the conditions for the maximum of the 
emission of a pair of photons, further work on the inter 
arrival times between such pairs and their entanglement 
is also of interest. 
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